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Abstract 

In baseball games, the coefficient of restitution of baseballs strongly affects the flying distance of 
batted balls, which determines the home-run probability. In Japan, the range of the coefficient of 
restitution of official baseballs has changed frequently over the past five years, causing the number 
of home runs to vary drastically. We analyzed data from Japanese baseball games played in 2014 
to investigate the statistical properties of pitched balls. In addition, we used the analysis results to 
develop a baseball-batting simulator for determining the home-run probability as a function of the 
coefficient of restitution. Our simulation results are explained by a simple theoretical argument. 
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I. INTRODUCTION 


The bounce characteristics of baseballs have a large influence in baseball games; thus, 
baseball organizations often establish rules concerning official balls. For example, in Major 
League Baseball (MLB), the baseballs are made by tightly winding yarn around a small core 
and covering it with two strips of white horsehide or cowhide^. 

For the estimation of the bounce characteristics, the coefficient of restitution e is widely 
used, which is dehned as 

"" f/F ^ ^ 

where Vi and Vj. are the speeds of incidence and rebound, respectively, in a head-on collision of 
a ball with a plane. Note that the coefficient of restitution determines the loss of translational 
energy during the collision. The coefficient of restitution depends on the kind of material 
and the internal structure of the ball, as well as other factors such as the impact speed^^^, 
impact angle^F^ temperature of the balF’®^ and humidity at which the balls are stored^°. 

Various baseball organizations officially determine the range of the coefficient of restitu¬ 
tion of baseballs. For example, the coefficient of restitution of an MLB baseball is required 
to be 0.546 ± 0.032^^. Regarding Japanese professional baseballs, the Nippon Professional 
Baseball Organization (NPB) hrst introduced their official baseball in 2011, which was used 
in both Pacific and Central League. Table I shows a chronological table indicating the av¬ 
erage coefficient of restitution for baseballs used in Japanese professional baseball games 
from 2010 to 2013^^, along with the annual number of home runs^^. Clearly, the number 
of home runs decreased drastically in 2011 compared with 2010, although the difference in 
the average coefficient of restitution is only on the order of 10“^. The average coefficient 
of restitution increased in 2013 because the NPB made a baseball equipment manufacturer 
change the specihcation of the baseballs in order to increase the level of offense in baseball 


games. 

Generally, the number of home runs is strongly affected not only by the coefficient of 
restitution of baseballs, but also other various factors, such as the climate, the specihcations 
of bats, and the batting skills of players, and so on. Sawicki et ah constructed a detailed 
batting model incorporating several factors, including the air resistance, friction between 
the bat and ball, wind velocity, and bat swing^^. Although they investigated the optimal 
strategy for achieving the maximal range of a batted ball, they did not calculate the home- 
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TABLE I. Coefficient of restitution of baseballs and number of home runs in Japanese professional 
baseball games. 


Year 

Coefficient of Restitution (average) 

Number of Home Runs 

2010 

0.418 

1,605 

2011 

0.408 

939 

2012 

0.408 

881 

2013 

0.416 

1,311 


run probability, because it may be difficult to choose proper parameters for the home run 
probability function. However, a quantitative research on the relationship between the 
coefficient of restitution of baseballs and the home-run probability is valuable for two reasons. 
First, Table. I indicates that the home-run probability strongly depends on the coefficient of 
restitution of baseballs because the small amount of changes in the coefficient of restitution 
can alter the flying distances of batted balls^®. Second, the coefficient of restitution of 
baseballs is a controllable factor that is important for the design of baseball equipment. The 
home run probability as a function of the coefficient of restitution can be a simple criterion 
to evaluate the characteristics of official baseballs. In addition, a quantitative research on 
the relationship between the coefficient of restitution of balls and the home-run probability 
is also valuable for physics education, as the problem is closely related to topics covered in 
undergraduate physics. 


In this study, we developed a batting simulator using real baseball data to quantitatively 
investigate the home-run probability as a function of the coefficient of restitution. This paper 
is structured as follows. In the next section, we describe the data analysis and analysis 
results. Sections 3 presents the construction of our batting simulator and the simulation 
results. In Sections 4 and 5, we discuss and snmmarize our results. Appendices A and B 
are devoted to the derivation of the averaged force in a binary collision between a ball and 
a bat and the algorithm for the collision, respectively. 
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II. DATA ANALYSIS 


To construct our batting simulator, we first analyzed pitching data for Japanese profes¬ 
sional baseball games held in 2014. We used data from Sportsnavi^®, which show various 
data about the pitched balls in an official game, including the ball speed, pitch type, and 
position of a ball crossing the home plate. Figure 1 shows a schematic of a part of a Sport- 
snavi page. In the data, the pitching zone is divided into 5x5 grids from the pitcher’s 
perspective, wherein 3x3 grids, represented by thick lines, corresponds to the strike zone 
(see the left panel of Fig. 1). The numbers and the symbols in a grid respectively show 
the order and types of pitches, respectively, at different positions on the grid. Information 
about each pitch, including the ball speed, is presented in the table shown in the right side 
of Fig. 1. For a later discussion, we numbered the horizontal and vertical positions of each 
grid as shown in Fig. 1. 

Using the Sportsnavi database, we manually recorded all the positions and ball speeds of 
pitches in 12 selected games held in Nagoya Dome Stadium in Nagoya, Japan, from August 
6 2014 to September 25 2014. We chose games held in indoor domes because the flight of 
baseballs is hardly affected by climatic factors such as the wind strength. We collected and 
analyzed data for 1,548 pitched balls. 
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1 2 3 4 5 

Horizontal index 


Number Type of Pitch Ball Speed 

1 curve ball 109 km/h 

2 straight fastball 141 km/h 

3 straight fastball 142 km/h 

4 two-seam fastball 143 km/h 


FIG. 1. Schematic of a part of the Sportsnavi page. 


Figure 2 shows the distribution of the pitched-ball speed n, where the open circles in¬ 
dicate the calculated probabilities as a function of v. To obtain the distribution function 
approximating these data, we divided the ball-speed data into two categories: those for 
straight balls and those for breaking balls having a curve, a two-seam fastball, etc. For each 
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of the categorized data points, we fit the normal distribution dehned by 

fiiv)= J- expj--^ — (2 = 1,2), (2) 

V27rcri I cr^ J 

where /ij and cxj are the mean and the standard deviation, respectively. The htting parame¬ 
ters are presented in Table II, where 2 = 1 and 2 = 2 correspond to the straight and breaking 
balls, respectively. 



FIG. 2. Distribution of ball speeds. Open circles show the probabilities at each ball speed. Solid 
black curve shows Eq. (3) with the fitting parameters shown in Table II. Solid red and blue curves 
show the distributions of the straight and the breaking balls, respectively, weighted with p = 0.45. 

Considering fi(v) (i = 1, 2) to be components, we hnally obtained the mixture distribution 
of the pitched-ball speed v as 

</>H=pfiM + (^-p)f2(v) (o<p<i). (3) 

Here, p is the mixing parameter. The black solid curve shown in Fig.2 indicates Eq. (3) 
with p = 0.45, which closely approximates the distribution of the pitched-ball speed with 
the coefficient of determination equal to 0.9942. The red and the blue curves show the hrst 
and the second terms, respectively, in the right-hand side of Eq. (3). Generally, the value of 
p represents the probability of selecting each component of the mixture distribution. Thus, 
we consider that the pitchers chose straight and breaking balls with probabilities of p = 0.45 
and 0.55, respectively. 

Next, we show the distributions of the horizontal position y and the vertical position 
of the pitched balls in the pitching zone. Figure 3 shows the distribution of the horizontal 
position of the pitched balls, where the horizontal axis indicates the positions of the grids 
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TABLE II. Fitting parameters used in Eq. (2). 


fJi [km/h] 

Pi [km/h] 

(72 [km/h] 

fi 2 [km/h] 

139.1 

4.58 

127.8 

7.80 


in Fig. 1, which are numbered from left to right. The open circles indicate the calculated 
probabilities as a function of the horizontal position. To obtain the distribution function 
approximating these data, we htted the mixture distribution of the normal distributions, 
which is expressed as 

= qhiiy) + {I - q)h2{y) (0 < g < 1), (4) 

where hi{y) and h 2 {y) are the normal distributions with the same standard deviation: 0.77. 

The means of hi{y) and h 2 {y) are y = 2 and |/ = 4, respectively. In Fig. 3, the black 
solid curve shows Eq. (4) with the mixing parameter q = 0.52, whereas the blue and red 
curves show the hrst and second terms, respectively, in the right-hand side of Eq. (4). The 
coefficient of determination of the ht is equal to 0.986, which indicate a good ht. This 
indicates that the pitchers tended to choose the right and left sides of the strike zone with 
a similar probability to prevent home runs. 



FIG. 3. Distribution of the horizontal position of pitched balls at the home base. Solid black curve 
shows Eq. (4) with q = 0.52. 

On the other hand. Fig. 4 shows the distribution of the vertical position for the pitched 
balls. The probability is almost constant, except at z = 5, which is the highest position. 
This indicates that the pitchers tended to choose heights all over the strike zone with a 
similar probability. The frequency at = 5 is relatively small to avoid the risk of long ball 


6 














hitting. On the other hand, the frequency at = 1 is almost same as those at z = 2,3, A 
inside the strike zone, which may be attributed to that the balls at z = 1 are difficult to 
strike. 



1 2 3 4 5 

Vertical Index 


FIG. 4. Distribution of the vertical position of pitched balls at the home base. 


III. SIMULATION 

We performed a batting simulation according to the results presented in the previous 
section. Figure 5(a) shows our simulation setup, in which a home plate was placed at the 
origin of a Cartesian coordinate system. The initial position of the center of mass of a 
baseball was set at tq = (18.44 m, 0 m, 1.8 m) according to the official baseball rules in 
Japan and the average height of Japanese professional baseball players. The baseball was 
pitched in the negative direction of the x axis with an initial velocity of Vq = VqC, where 
C is the unit direction vector of the pitch, which will be defined later. Here, Vq is randomly 
chosen according to the distribution function given by Eq. (3) with the parameters shown 
in Table II. 

The pitch direction was defined as follows. A ball needs to be thrown from about shoulder 
height and be launched in an almost horizontal direction to cross the plate at the correct 
height^. Thus, we assume that a pitcher throws a ball toward a point P = (0, Py [m], 1.8 m) 
that is on the y-z plane. A thrown ball travels a curved path to cross the plate at a height 
less than the initial height due to the gravitational force (Fig. 5(b)). Py is the random 
variable selected by the distribution function ipiPy) defined in Eq. (4). The unit vector of 
the pitch direction was defined as C = C/|C|, where C = P — fq. 
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FIG. 5. Schematics showing (a) our simulation setup and (b) the definition of 

A. Modeling of Bat and Swing 

We consider the bat to be an uniform cylinder 1 m long and the diameter of the base to 
be 6.6 X 10“^ m, in accordance with the official baseball rules. The bat is placed along the 
y axis, and its center of mass is positioned at (0, 0, Zb [m]), where Zb is defined later. For 
later discussions, we respectively label bases of the bat as A and B so that the y component 
of the center of B is larger than that of A. 

Let us assume that a thrown ball crosses the y-z plane between the time t and t + At, 
where At = 0.05[s] is the time step of our simulation (Fig. 5(b)). First, we calculate using 
the X component of the intersection between the y-z plane and the line segment connecting 
r(t) and r(t + At), where r(t) is the position of the center of mass of a thrown ball at the 
time t. Next, we determine Zb as follows: 


Zb = z[ + (Jb, (5) 

where ab represents the random numbers chosen from the normal distribution with a mean 
of 0 and a standard deviation of 0.0366 m, which is the diameter of an official baseball in 
Japan. 


CO 



A 


FIG. 6. A schematic of a collision between a ball and a bat. 














Figure 6 shows a schematic of a collision between a ball and a bat from the viewpoint of 
the home base. We assume the bat swings around an axis passing through the center of B. 
The bat swings with the angular velocity lJ = (—sin6', 0, cos 6*) with = 34 rad/s, 
which is the typical value of bat swings. Here 9 is the angle between the direction of uj' and 
(jj = (0, 0, Uz), the value of which is randomly chosen from the range 0° < 0 < 10°. When a 
ball collides with the bat, we calculate the vector f, which is dehned by the vector from the 
center of B to the foot of the perpendicular line passing through the center of the ball to the 
central axis of the bat. The velocity of the bat on the collision is defined by V;, = a;' x f. 


B. Equation of Motion of Ball 


Basically, a pitched ball obeys the equation of motion 

CpT 

+ Fl + mg. 


( 6 ) 


where m is the mass of the ball. The three terms on the right-hand side of Eq. (6) represent 
the drag force, the magnus force, and the gravitational force, respectively. 

The drag force F^) generated by the resistance from the air decreases the ball speed. F 
is expressed as 

Fd = ^CdpV^AY, (7) 

where Cd, p. A, and V are the drag coefficient, the density of air, the cross section of a 
baseball, and the unit vector in the direction of the velocity, respectively. We use Cd = 0.4, 
corresponding to a baseball pitched at a high speed of approximately 40.2 to 44.7 m/s^. Our 
choice of the value will be discussed in later section. In addition, we use p = 1.29 kg/m^, 
which is the value at 0 C° and 1 atm. 

On the other hand, the magnus force F^ generated by the backspin of the baseball 
enhances the flight of the ball after it is thrown and batted^^’^®. F^ is expressed as 

Fl = ^CLpV^AVn, (8) 

where Cl and are the lift coefficient and the unit vector perpendicular to V, respectively. 
We used Cl = 0.2, which is a typical value for a spinning balT. The gravitational acceleration 
is indicated by the vector g = (0,0, —9.8 m/s^). We numerically solved Eq. (6) using the 
velocity Verlet scheme^®. 
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C. Condition for Collision 


A thrown ball changes its flight direction when it collides with a bat. We assume that 
the ball collides with the bat when all of the following conditions are fulhlled: 

1. The distance between the center of mass of the ball and the central axis of the bat is 
less than the sum of the radius of the ball and that of the bat. 


2. The y component of the ball’s position, r^, satishes —0.5 m < < 0.5 m. 


When these conditions are fulhlled, the ball is rehected by the averaged repulsive force 
during At, 


F 


yil + e)V ■ n 

-1-n 

At 


( 9 ) 


where /r, n and V are the reduced mass of the ball and the bat, the normal unit vector that 
is perpendicular to the tangential plane between a colliding ball and a bat, and the relative 
velocity of the ball to the bat, V = V — V^, respectively. 

In Eq. (9), we use the coefficient of restitution between a ball and a bat e, which is 
dehned by 


e = e 



m 

Afeff 


+ 


m 


( 10 ) 


where Mgff = I jr with the moment of inertia I of the bat around the center of The 

derivation of Eq. (9) and the algorithm for the collision are summarized in Appendices A 
and B, respectively. 

When a batted ball falls on the ground, we calculate the distance D between the point 
of landing and the home plate. 


IV. SIMULATION RESULTS 

Figure 7 shows the probability densities of the flying distance D of batted balls (e = 0.41 
and e = 0.45), each of which was calculated from 1,000 samples. The highest peak position 
shifts from 75 m to 125 m with increasing coefficient of restitution of a ball. The total 
frequencies in D > 150, which can be regarded as the number of home runs, increases 
with increasing coefficient of restitution indicating an increase in the home-run probability. 
Thus, we investigate the relationship between the home-run probability and e by changing 
the value of e from 0.41 to 0.45. 
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FIG. 7. Probability density of flying distance of batted balls. Solid and dashed curves show the 
results of e = 0.41 and e = 0.45, respectively. 



FIG. 8. Relationship between e and home-run probability P(e}. Solid curve shows Eq. (13) with 
Cl = 242 m and a = 52.2 m. 


Figure 8 shows the relationship between the home-run probability and e. Here, each 
data point was calculated from 1,000 samples. We consider a sample with D > 150 m as a 
home run and dehne the home-run probability as the ratio of the number of home runs to 
1,000. As shown in Fig. 8, the home-run probability increases with e, which is intuitively 
understood. 

Here, we estimate the functional form of P{e) using a simple theoretical argument. Let 
us suppose that a projectile is launched from the ground with a launch speed v' and a launch 
angle do under the air friction. The range of the trajectory under the air friction can be 
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estimated by a linear function of v (Ref. 2), so that we roughly estimate the range of the 
trajectory as a linear function of the coefficient of restitution e as Cie with the constant Ci. 
Here we assume that the speed of the thrown ball upon the collision with the bat is almost 
constant. 

We assume that the probability density of the flying distance of a batted ball which will 
land at approximately D = 150 m can be approximated by the normal distribution as 


p{D) = 




: exp 


{D - 

2a2 


By integrating Eq. (11) from D = 150 m to inhnity, we obtain the probability function 


P(e) = / p(P)dR) 


'150 


1 — erf 


150 - Cxe 


V2i 


a 


( 11 ) 

1 

( 12 ) 

(13) 


where erf(a;) is the error function. We assume both Ci and a as htting parameters. The 
solid curve in Fig. 8 corresponds to Eq. (13) drawn with a = 59.8 m and Ci = 271 m, 
which indicates a good estimation. 


V. DISCUSSION 

Figure 3 shows the distribution of the horizontal index y of pitched balls at the pitching 
zone, which can be approximated by the mixture distribution of the normal distributions. 
Because of the method of division of the pitching zone in the database that we used, the 
number of points was too scarce to identify a more detailed distribution. To obtain a more 
detailed distribution of the position of pitched balls, we can use more accurate data from 
PITCHf/x^. However, from the viewpoint of pitcher strategy, it may be desirable to aim at 
the edges of the strike zone; thus, we may hud two peaks at the right and the left edges of 
the strike zone in the real distribution. 

Figure 7 shows the probability density distributions p{D) of the distance D, each of 
which has some peaks. For comparison with real data, a Japanese group attempted to draw 
a histogram of the distance of batted balls using data from Ultimate Zone Rating (UZR)^^ in 
2009 and 2010. Their results show that the histogram has two peaks around 40 and 100 m, 
which is qualitatively close to our simulation results at e = 0.45. However, their results show 
that the probability of Ending a ball at the distance 40 < H < 80 is very low, although our 
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simulation results show that the probability density at e = 0.45 is almost same in the same 
range. Notice that the group constructed the histogram using the positions at which the 
helders caught or picked up the batted balls, which are sometimes different from the points 
of landing. Thus, their histogram can be considered as the mixture of the histograms of the 
position of pickup by the inhelders and the outfielders, which may cause the discrepancy 
between their results and our results. 

In the UZR analysis, the probability of hnding a batted ball between 91.44 and 121.92 m 
was 0.179. In our simulation, on the other hand, the home-run probability P{e) was ~ 0.17 
at e = 0.42, which is the closest value to the averaged coefficient of restitution in 2010. The 
discrepancy between the two values is attributed to that the data in 2009 and 2010 were 
included in the UZR analysis. 

Finally, our model is a simplihed model based on simple rigid-body collisions, which 
ignores the rotation of a batted ball. Our model includes the lift coefficient Cl = 0.2 as a 
hxed value, although Cl can vary depending on the rotation of the balh’^^. For example, 
recent experiments have demonstrated that Cl strongly depends on the spin parameter, S = 
Ru/V, where u and R are the angular velocity and the radius of a ball, respectively^®’^^’^^. 
On the other hand, Cd also varies according to Reynolds number, although our model 
includes the drag coefficient Cd = 0.4 as a hxed value. Especially, the drag force on a 
baseball is known to drop sharply at speeds typical for thrown and batted balls^®’^^’^^’^^. 
To improve our model, we will need to incorporate Cl and Co dependant on the motion 
of thrown and batted balls. In addition, constructing a model based on elastic collisions in 
which the friction between the ball and bat is considered will yield more accurate results. 


VI. CONCLUDING REMARKS 

In this paper, we analyzed real data for the speed and the course of pitched balls in 
professional baseball games in Japan. Our results show that the distribution of the ball 
speed can be approximated by the mixture distribution of two normal distributions. In 
addition, we found that the horizontal position of pitched balls in the pitching zone obeys 
the mixture distribution of two normal distributions, each of which has a peak at the edge 
of the strike zone. 

We simulated collisions between baseballs and a bat, where the statistics of the pitching 
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obey our analyzed results. We finally obtained the probability density distribution p{D) 
of the distance D between the home plate and the landing point of the batted balls to 
calculate the home-run probability as a function of the coefficient of restitution. By using 
a simple theoretical argument with the assumption that p{D) around D = 150 m can 
be approximated by the normal distribution, we quantified the home-run probability as a 
function of the coefficient of restitution e of baseballs. 

As stated in the previous section, we will obtain more accurate results in future works 
by improving our model. Developing a methodology to calculate the home-run probability 
may yield useful information for designing baseballs. 

Appendix A: AVERAGED FORCE IN BINARY COLLISION OF RIGID BODIES 

In Ref.l, the author derived the rebound velocity of a ball in a head-on collision with a 
bat^. In this appendix, we derive the rebound velocity of a batted ball in a three dimensional 
oblique collision. Note that we ignore the rotation of a ball and a bat in our argument. 

Figure 9 shows a ball of mass m colliding with a bat of mass M. Here we denote the 
colliding velocities of the ball and the bat as V and V^, respectively. Assuming that the 
ball experiences the averaged force F from the bat during the duration At, we can describe 
the rebound velocities of the ball and the bat, V' and as 

V' = V + —At, (Al) 

m 

v; = Vi - (A2) 

from the definition of the averaged force and Newton’s third law of motion. Note that the 
total momentum of the system is conserved after collision, MV;, -|- mV = MV^ -|- rnV'. 

By subtracting Eq. (A2) from Eq. (Al) and introducing the reduced mass p defined by 
l//i = l/m-|- 1/M, we obtain 

V' = V + -FAt, (A3) 

/i 

where V and V' are the relative velocities of the ball to the velocity of bat before and 
after collision, respectively. By introducing the normal unit vector n perpendicular to the 
tangential plane of the bat and ball, the scalar projection of V' onto n is calculated as 

V' • n = V ■ n + -FAt, (A4) 
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FIG. 9. A schematic of a collision between a bat of mass M and a ball of mass m. The velocities 
of the centers of mass of the ball and bat are denoted by V and V;,, respectively, n is the normal 
unit vector perpendicular to the tangential plane. 

where we used F = Fn. 

Using Eq. (A4) in the definition of the coefficient of restitution e between ball and bat, 

|V-n| =e|V-n|, (A5) 


we obtain the equation 


/3 +V-n(l 



/3 + V ■ n(l + e) 


0 . 


(A6) 


where {3 = FAt/fi. The two solutions of the quadratic equation Eq. (A6) are respectively 
written as 


(3i = -V ■ n(l - e), 
132 = -V -11(1 + 6). 


(A7) 

(A8) 


In the two solutions, only /32 corresponds to the averaged force in the collision because I3i 
becomes 0 (F = 0) when e = 1, which will cause the penetration of the ball into the bat. 
Thus, the averaged force is calculated as 


F = Fn 


h /^2 


-n 


At 

+(1 + e)V ■ n 

aF 


-n. 


(A9) 

(AlO) 
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FIG. 10. A schematic of a collision between a bat (large circle) and a ball (small circle). The ball 
at time t will penetrate into the bat at the time t + At, which is unrealistic. To obtain the normal 
force acting on the ball, we need to calculate the time t + At when the ball touches the bat. 

Appendix B: ALGORITHM FOR COLLISION BETWEEN BALL AND BAT 

In appendix A, we explained the way to calculate the averaged force acting on the ball 
during collision. However, it is difficult to obtain the unit normal vector n in naive calculation 
because the ball can penetrate into the bat during the simulation time step, which is due to 
the nature of the hnite difference approximation of derivatives. In this appendix, we explain 
the algorithm of the collision between a bat and a ball to avoid penetration used in our 
simulation. 

Figure 10 shows a schematic of a collision between a ball and a bat, where a ball at the 
time t is going to collide with the bat. When the hrst condition of collision (see section 
III.C) is fulhlled, the ball has penetrated into the bat at the time t + At. Thus, we put the 
ball back to the previous position at the time t to determine the remaining time At' before 
collision. From the condition that the distance between the center of the bat and the center 
of the ball at the time t + At' equals to the sum of the radius of the bat and that of the ball, 
we obtain a quadratic equation for At'. Here we assume that the ball travels linearly with 
a constant velocity during At' to simplify the calculation. Among the two solutions for the 
quadratic equation, we choose the positive one for At' which has a physical meaning. 

From the position of the ball at the time t + At', we can calculate all the variables such as 
n, f, Vf, and e in Eq. (9). The ball reflects from the bat according to the following algorithm. 
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First, we put the ball back to the position at the time t. Next, we apply the half of the 
averaged force to the ball during 2At so that the impulse from the bat becomes constant. 
With this two-step time evolntion, the ball bounces from the bat without penetration. 
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